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PACS 03 . 65 . Pm - Relativistic wave equations 

PACS 03.65.Vf - Phases: geometric; dynamic or topological 

PACS 03 . 65 . Sq - Semiclassical theories and applications 

Abstract - The recent discovery of electron vortices opens up a wide research domain previously 
unexplored. The present paper explores the relativistic properties of these singular fermion beams, 
and quantifies deviations from classical wave theory. It is common in electron optics to use the 
Schrodinger equation neglecting spin. The present paper investigates the role of spin and the 
total angular momentum Jz and how it pertains to the vortex states. As an application, we also 
investigate if it is possible to use holographic reconstruction to create novel total angular momen- 
tum eigenstates with a space- dependent spin polarization in a Transmission Electron Microscopy 
(TEM) setup. It is demonstrated that relativistic effects disappear in the paraxial limit, and spin 
effects in holographically created electron vortex beams can be safely ruled out. 
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Introduction. — Singular optics is the study of light 
waves with a phase singularity ^LJ. For two decades, since 
the experimental discovery of optical vortices [2 , study of 
singular properties in the scalar theory of light has led to 
fascinating applications in various areas of physics [3]-[8]. 
Much more recently, the concepts and experimental se- 
. tups from singular optics found their way to electron mi- 
' croscopes [QplOj . Different methods are being investigated 
for the production of electron vortices pT|[T2] and possible 
applications at the nanoscale are surfacing [13]. The first 
steps into developing a theory of electron vortices have 
been taken p!4Hl6] The quantum mechanical description 
. of electron optics coincides for the most part with scalar 
optics. This is not surprising due to the fact that for sta- 
tionary states, both the Maxwell and Schrodinger equa- 
tions take on the exact same form. 

Both photon polarization and electron spin determine a 
spin angular momentum (SAM) vector, which contributes 
together with the orbital angular momentum (0AM) vec- 
tor to the total angular momentum (TAM) vector: J = 
L -\- S. In this paper, the cylindrically symmetric case is 
studied, and the quantities of interest are the projections 
of the angular momenta on the z-axis (parallel to beam 
propagation): Jz = Lz -\- Sz- In what follows, cylindri- 
cal symmetry is assumed, and all references to "angular 
momentum" refer to the projected angular momenta on 
the vertical z-axis. One can pose the question of the in- 
fluence of the bosonic nature of the photon versus the 
fermionic character of the electron. Electrons have half- 



integer spin, which translates to half-integer TAM as a 
proper observable, where photons have integer SAM, and 
thus integer TAM. The difference in TAM is elucidated 
by inspecting the Dirac equation, for which Lz is not a 
good quantum number, and it is only the TAM Jz that 
commutes with the Dirac Hamiltonian. It may therefore 
come as a surprise that the scalar theory works so well for 
fast electron waves. Yet experiments show that electron 
vortex beams behave as their scalar counterparts, and e.g. 
spin seems to have little to no influence in the entire TEM 
(Transmission Electron Microscopy) image formation pro- 
cess pTHTO] . depending on the acceleration voltage used. 
This is caused by the relatively small magnitude of the 
spin effects (or conversely, the high electromagnetic fields 
required to show the effects), and here it is shown that in 
the paraxial limit one can speak of proper Lz eigenstates 
even in Dirac theory. 

Once electron spin comes into play, there are two scalar 
components (at least in the non-relativistic Pauli case) 
that cooperatively determine the overall quantifiable prop- 
erties of the beam. These components are independent if 
there is no external magnetic field. This leads to the pos- 
sibility of distinct Jz eigenstates, which possess interesting 
properties such as a skyrmionic spin distribution [20J . We 
pose the question whether such Jz states can be produced 
in similar ways as the Lz eigenstates, e.g. using holo- 
graphic masks [HITO]. In what follows, natural units are 
used: h = c = 1. Operator expectation values are de- 
noted (O). When the state is an eigenstate, its values are 
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denoted without brackets as O. 

Phase singularities and electron vortices. — The 

simplest description (using a scalar wave theory) of a phase 
singularity in the paraxial regime in a transverse plane can 
be written down in terms of the wave function (or in optics, 
the electric field E)^ with a singular phase factor: 



^V^(r). 



(1) 



One can readily see that this is an 0AM eigenstate: 
= —id^p^ = m^. The phase increases linearly with 
ip^ resulting in a singular point at r = in the center 
of the circular phase ramp, where the phase is undefined. 
Experimental production of these types of wave functions 
can be accomplished by using computer-generated holo- 
grams . These binary masks can be calculated from the 
interference of a tilted plane wave with the target vortex 
mode. Using the principles of holographic reconstruction, 
illuminating the mask with the reference beam generates 
the target wave as output. In the case of binary masks, 
one also gets several higher diffraction orders in the far 
field. This scalar optics technique allows to create any 
0AM and even superpositions of Lz are possible. Exactly 
the same technique was used with fast electron waves cre- 
ating specific 0AM states and superpositions in ^21]. To 
understand why a scalar theory apparently is sufficient 
for describing particles with spin we will first consider the 
non-relativistic Pauli equation. 

OAM Pauli spinors. — 

Mathematical solutions. The non-relativistic equa- 
tion of choice to describe electron polarization is the 2- 
component Pauli equation. It provides a non-relativistic 
description of spin and its interaction with electromag- 
netic fields, which describes the processes in an electron 
microscope sufficiently. The (stationary) Pauli equation 
(assuming ip oc ^f{r)e~^^^) with an electromagnetic field 
can be written as follows, 

1 

2m 



cr{p- eA)Y + + (Ibct - B]^! = E^, (2) 



with vector potential A = (^,A) representing an exter- 
nal electromagnetic field, S = V x A, the scalar (elec- 
trostatic) potential, a the Pauli vector, e the electron's 

charge, I-^b = ^ the Bohr magneton, and ^ = ( ^2 j the 

2-component Pauli spinor. In cylindrical symmetry ex- 
pressed in the coordinates (p^cp^z)^ the general solutions 
are readily given by 







(3) 



These states are eigenstates of the energy E, the forward 
momentum /c^, the transverse momentum the OAM 



Lz = n, and the spin Sz (labeled with the ± superscript), 
and thus also of Jz = Lz -\- Sz- The radial functions Jn{x) 
are the n-th order cylindrical Bessel functions of the first 
kind. In principle, one can combine any ^+ with any other 
and still have a valid solution, although in general it 
won't be an Sz or Lz eigenstate anymore. The case where 
m = n + 1 is interesting because this combination is also 
an eigenstate of the TAM Jz = n 
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±e^(-+i)^Jn+i(/^p) 



(4) 



This state has traded Lz and Sz for Jz and the non- 
relativistic limit of the transverse helicity h± = a ■ ^^f^- 
These states are not OAM eigenstates, and the two scalar 
spin components have a different OAM. This state's Lz 
expectation value is equal to its Jz eigenvalue: 

(Lz) = ^[n + (n + 1)] = n + i = J„ (5) 

This should come as no surprise because {Jz) = (Lz) + 
(Sz) and (Sz) is equal to for a spin unpolarized beam. 
Since the two-component Pauli theory poses no helicity 
constraints, the transverse spin vector orientation can be 
defined by choosing an appropriate phase angle between 
the spin up and spin down component. This can result 
in a net spin polarization in the transverse plane directly 
resulting from this relative phase. 

Selection on Lz by a holographic mask. Fraunhofer 
diffraction of a scalar function has been studied extensively 
in the optics literature using the stationary wave equation: 



V'^E^k^E = 0. 



(6) 



This Helmholtz equation can be derived from the quan- 
tum mechanical Schrodinger equation when a stationary 
state is considered, leading to the same Kirchhoff diffrac- 
tion formalism applying for light and matter waves. In the 
case of the 2-component Pauli equation, one has to inter- 
pret its spin-less limit leading to the Schrodinger equation. 
This has been considered formally previously in [22 , and 
the Pauli equation is found to reduce to the Schrodinger 
equation with the proper probability current only if one 
considers a spin eigenstate without its interactions with 
a magnetic field. The effective decoupling of both Pauli 
spinor components in the absence of an electromagnetic 
field into two non-interfering components, is entirely con- 
sistent in this case and can be described by the following 
relation: 



V 



5(V^^)0 5(V^') 



(7) 



where V signifies the Pauh and S the Schrodinger equa- 
tion. A spin unpolarized beam diffracting on an Lz binary 
fork mask can be described by two seperate scalar wave 
optics systems, each selecting on the same OAM. One ef- 
fectively has two independent and non-interfering wave 
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functions undergoing Fraunhofer diffraction, leading to an 
unpolarized vortex state ^ = (^+ + for unpolar- 
ized incidence (see ([3])). 

A general Pauli electron vortex beam can be parametrized 
with a G R (the radial functions /(r) and g{r) are as- 
sumed to be normalized to 1): 



^(r, if) 



(8) 



where n and n' are the two spin components' 0AM. This 
normalized state (|^P = 1) is an exact solution of the 
field-free Pauli equation, but not (in general) an eigenstate 
of any angular momentum. 

If the weight of both components is equal (a = 1), one has 



{Lz) = ^(n + n), 
which leads to an = n eigenstate if n = n', and 



(9) 



(J,) = \[{n + 1/2) + (n' - 1/2)] = \{n + n'), (10) 

which leads to (J^) = n if n = n^ So in the case 
of identical vortices in both spin components, we have 
Lz = {Jz) = ^- In terms of measurable effects, the system 
behaves as two independent and identically propagating 
scalar systems in an Lz eigenstate. 

If the weight of both components is not equal, one has for 
the angular momenta expectation values: 



{Lz) = 



n - 
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- 1, spinor ([8]) is an eigenstate 
^, and this is unaffected by 



In the special case n' = n 
of Jz with eigenvalue n H 
the value of a. The above shows that even though only 
Jz is a proper quantum number, for the special case of 
unpolarized electrons, the Jz eigenstate also has {Lz) = 

Jz- 

Fig. [T] shows the scalar far field image of an illumi- 
nated fork mask, showing clearly the 0th and first order 
vortices (±1 and ±2). The third far field image shows the 
density of a spinor Jz eigenstate (j4|), which is the sum of 
the other two images. It contains two superimposed and 
independent scalar vortices of different order. 

Selection on Jz by a holographic mask. In order to 
create a Jz eigenstate such as those in (|4]), one needs a 
superposition of different 0AM for each spin component. 
For photons, subwavelength gratings can be used to pro- 
duce such states [23]. This is difficult for electrons in an 
electron microscope, as the wavelength of an accelerated 
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Fig. 1: The farfield of a scalar n = 1, n = 2, and a spinor 
Jz = ^ (equation Q) with the radial distribution caused by 
a circular aperture with a fork pattern. The doughnut shapes 
are phase vortices, with a phase ramp of to n27r, and a cen- 
tral zero-density phase singularity. Intensity is saturated which 
allows direct comparison with fig. [2] 



electron is typically of the order of pm. Creating half- 
integer holographic masks does not have the desired ef- 
fect. As previously shown by Berry [24], this leads to 
non-integral 0AM which is nothing but the sum of inte- 
gral 0AM contributions, which show in the wave pattern 
as the interference of the different 0AM modes. A holo- 
graphic mask designed for half integer Lz will diffract both 
spin components in exactly the same way, leading to no 
space- variant polarization as described by Wang et al. [20] 
for the target states. The relative phase between the two 
spin components is important for the interference though, 
as shown in fig. [2] and explained below. 
To demonstrate this, we calculated a binary holographic 
mask from the interference pattern of an unpolarized (in 
z-direction) incoming plane wave, and a state of the form 
(|4|). Fig. [2] shows several possible masks for m = n + 1 
and different relative phases (=transverse spin orienta- 
tion) between the two spin components of the incoming 
tilted wave. The reference tilted plane wave used is: 



Ref - 



(12) 



where C is equal to 1, — z, and ^{1 — i) for the three figures 
in fig. [2j As one can clearly see in fig. [2l the relative phase 
factor between the two spin states (which determines the 
direction of the transverse spin vector) used to calculate 
the holographic mask, also determines the orientation of 
the additional dislocation line and the direction in which 
the extra integral order vortices split off from the most 
intense vortex. Once the mask is calculated, it does not 
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relative phase: relative phase 



relative phase: 




Fig. 2: Above: Holographic binary masks calculated by in- 
terfering Jz eigenstates with different relative phases for the 
incoming tilted spin states ([T2j). One can see that the relative 
phase influences the orientation of an additional dislocation 
line within the mask. Below: Far field images obtained by illu- 
minating the holographic masks with the tilted reference wave. 
The intensity is saturated to bring forward the structure of the 
first order vortices. Scalar diffraction gives the same result, 
and it is not equal to the third image in fig. [1] which shows 
the wanted Jz eigenstate. 



Equation 
the form: 



T3]) admits cylindrically symmetric solutions of 



V- 



-isk± 
[sk_\_ 



(14) 



E+ra 



These are eigenstates of kzj^ Jz^ and the transverse 

^ , which takes the values 



5^,3 ^ 



hehcity defined by h± — 7^7 

±1. Taking a similar linear combination as the one that 
relates the Pauli spinors in equations (j3j) and (j4]), one can 
construct from (p!4|) only an approximate eigenstate: 



^(+) = i 



n,s=+l ~t~ ^n,s= — 1) 



OC 







E+ra 



(15a) 



matter what the phase factor of the illuminating reference 
wave is; as shown above both components scatter indepen- 
dently, and the incoming beam's spin is of no influence to 
the resulting interference pattern. Due to considerations 
expressed in the previous section, this hybrid 0AM mask 
does not lead to the Jz vortex eigenstates. The result- 
ing wave is not of the desired form nor character, instead 
breaking down to a combination of integer order vortices. 
This shows the Jz eigenstates can only be obtained when 
we can influence both spin components independently (e.g. 
as with an electromagnetic field). Fourier reconstruction 
does give the expected result for spinors. 
We will now investigate under what conditions the intrin- 
sic relativistic spin coupling (detailed in the next section) 
might be exploited in the kinematical regime of a TEM. 



OAM Dirac spinors. — 

Mathematical solutions. In order to further under- 
stand the role of spin and thus the TAM of electron vor- 
tex states, the solutions of the cylindrical Dirac equation 
are considered. The underlying foundation of the ques- 
tion whether Lz or Jz is the most important variable, can 
be answered quantitatively by inspecting the fully rela- 
tivistically correct cylindrical solutions to the free Dirac 
equation (7^ are the Dirac matrices, ^ is a 4-component 
Dirac spinor): 

(z7^a^ - m) ^ = 0. (13) 



^i-l = \ (^n-l,.=+l - ^n-l,.=-l) 

e-^Jn(/c^r) 

oc _^,.(n-l).j_^(^^,) . (15b) 

V ^e^"^^n(/c^r) / 

The solutions in equation (fT5)) bear the closest resem- 
blance to the Pauli spinors discussed above (see (j3j)). 
These solutions are still Jz eigenstates with eigenvalue 
n ± 1/2, and become Lz = n and eigenstates only 
in the paraxial approximation: k±_ <^ E -\- m. Note that 
both of these solutions carry a different unapproximated 
eigenvalue of J^, although their approximate Lz eigen- 
value is equal. It is as though the paraxial limit removes 
the coupHng of the electronic wave function to its spin in 
this field-free case. Bliokh et al. [25] have previously pre- 
sented an alternative form of ([15]) , but the form here bet- 
ter captures the kinematical quantities that are relevant 
in a practical TEM experiment. The quantitative analysis 
of the relativistic spin coupling these solutions contain is 
presented below. 

Relativistic properties of electron vortex beams. As 
previously illuminated by Bliokh et al. there are two 
special states which, in principle, permit us to directly ob- 
serve a relativistic effect. There, it is shown that for the 
states and ^ (see ([T5|) ), there is a significant cen- 
tral contribution to the density p. Below is shown however 
that this effect is extremely small for typical TEM param- 
eters. The zeroth order Bessel functions are the only ones 
that are nonzero at the origin (i.e. the center of the beam). 
This gives the states ^^'^ and ^ (see equation ([T5|) ) a 
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Fig. 3: Central density distribution ([T6j) using the parameters 
mentioned in the text. The blue (solid) line is given by (|16p 
and the red (dashed) line is only the first term of that equation. 



nonzero central density: 
p(+^\r) = p[-\r) = 



(£; + m)2 



(£; + m)2 



Jo{k±r). 



(16) 



The third term results in a non-zero density at r = 0, 
proportional to - Approximating the radius of the 

vortex by only looking at the maximum of J^, one can 
estimate k± (up to O {{k±r)^)) for a single-/c^ state (see 
Appendix) : 



k± (in keV) 



0.37... 

R (in nm) 



(17) 



In state of the art TEM experiments, one can achieve a 
focused electron vortex [1^ of R ^ 0.5 A, and one has 
k± ^ 7.4 keV. For a 200 keV beam, this results in a central 
contribution to the density of about 3.7x10"^. This figure 
needs to be compared with the unity sized contribution 
of the other terms in (fT6|) . Fig. [3] shows a comparison 
of the first term of (p!6|) with the full expression for the 
aforementioned values. It should be compared with fig. 
2 in [25]. The relativistic effect for these states in this 
parameter regime are much smaller than the differences 
shown in that article. By increasing the focus and thus by 
further shrinking the vortex, one can in principle increase 
the magnitude of the effect, because k± increases with 
decreasing vortex size according to (pT|) . Though as one 



increases k±^ the relevant area of the central density will 
shrink as well. One can easily estimate the radius inside 
which the zeroth order contribution is larger than the first 
order in ([T6|) (see Appendix): 



rc 



0.24... 



^/c^ + m(m + Y^/c^ + m^^ 



1/4 



^/R 



(18) 



where R is in nm and the kinematic quantities in keV (nat- 
ural units). For our values of the parameters k_\_ and kz, 
one has rc ^ 4 pm, which is of the order of the electron 
wavelength. The area inside the critical radius is thus of 
the order of 50 pm^, in which the intensity is equal to ap- 
proximately 10 ~^ the intensity at the vortex maximum. 
For a spin-unpolarized beam the oppositely polarized vor- 
tex state (see (p!5|) ) is also present, which adds another 
unwanted unity sized background contribution which has 
to be discriminated upon measurement. The extremely 
small size of these effects allows us to say with certainty 
that pure relativistic effect of free-space electron vortices 
is unobservable for the conditions prevalent in an electron 
microscope. The smallness of these effects will translate 
into the fact that the Pauli description suffices, that is 
to say, ([7j) is valid, and effectively Lz can be maintained 
as a good quantum number even though this is not the 
case in the Dirac description. The extreme small size of 
the spin coupling effects in electron vortices lead us to 
conclude that single-electron diffraction through a holo- 
graphic mask cannot produce Jz vortex eigenstates. The 
Pauli equation supplies enough information for the TEM 
system and it results in per-spin Fraunhofer diffraction, 
leading to the coherent sum of the two spin states being 
present in the resulting diffraction pattern. 

Conclusions. — This paper has investigated the prop- 
erties of electron vortex states as they appear in electron 
microscopes. The question of the right quantum num- 
bers to qualify electron vortex states in a field-free region 
of space is answered. It is apparent that Lz is the right 
quantum number to label free electron vortex states, and 
neglecting spin in the description of their formation as long 
as no magnetic effects occur. This follows from the decou- 
pling of the non-relativistic Pauli equation in the absence 
of an electromagnetic field. In an attempt to quantify the 
field-free relativistic effects as previously illuminated by 
Bliokh et al. [25], an alternative form of the Dirac spinors 
representing approximate Lz electron vortex eigenstates is 
presented and their connection to the transverse helicity 
eigenstates of the Dirac equation is shown. These were 
used to calculate the size of the contribution of the rela- 
tivistic description to the central density present in specific 
states. It was found that for realistic conditions in modern 
electron microscopes, the field-free relativistic effects are 
very difficult, if not impossible, to observe: the relative in- 
tensity of this central density compared to the maximum 
of the full wave lies around 10 and the size of the area 
in which the effect dominates is of the order of the electron 
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wavelength. The current experimental apparatus does not 
lead to sufficiently large relativistic spin effects, so these 
cannot be measured. 

The usage of a holographic mask to poduce eigenstates 
was investigated. A holographic fork aperture works spa- 
tially, whereas spin is an additional internal degree of free- 
dom not directly related to the spatial diffraction phe- 
nomenon. The introduction of an electromagnetic field 
will couple the two components, and we expect the Jz 
eigenstates presented here to become important. One can 
thus not fabricate a Jz holographic mask without some 
form of magnetism, which could in itself lead to trans- 
verse spin-splitting as in [26], although it is not clear how 
the transverse Stern-Gerlach effect calculated there will 
couple to vortex formation by a holographic mask. A 
combined magnetized and forked pattern seems like the 
most achievable solution, once the necessary magnetiza- 
tion conditions can be experimentally achieved. Without 
magnetization or ultra-relativistic effects, one cannot in- 
fluence the spin degree of freedom by simple diffraction. In 
principle, one can repeat the calculations in this paper for 
neutrons, where perhaps the different interactions and pa- 
rameters can cause the effects described above to become 
observable. Neutron vortex states have not been presented 
as far as the authors are aware, but might bring the rich 
properties of phase singular beams to neutron optics. 
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Appendix. — 

Relation between the radial size of a Bessel beam in nm 
and k±_ in keV. In what follows no units are specified, 
because it is assumed has the inverse unit of r, resulting 
in a dimensionless argument to the Bessel function. An 
extremum of the first order Bessel function of the first kind 
is determined by the following equation: 



drJi{k±r) = 0, 



(19) 



which can be rewritten using the Bessel function's recur- 
sion relations to read: 



Mk±r) - J2ik±r) = 0. 



(20) 



The next step is expanding the Bessel functions up to 
order (/c^r)^ around k±r = (which turns out to be a 
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Fig. 4: Approximation of the radius of the vortex in function of 
k± for a Bessel beam of order 1. The curves are almost exactly 
superposed 

vahd approximation for k± > 1). This leads to following 
even fourth degree equation: 

(fc^ (fc^_(fc^_(fc^ 

4 + 64 " 8 96 ■ ^ ^ 

Solving this for t = ^^^^'> , gives the relevant root: 

^ = Vo:8835 ^r=if^. (22) 

2 k± 

The result is shown in figure HJ The maximum relative 
error in the domain k± G [1,10] of this approximation is 
under the 0.4%. 

The final step is the conversion to natural units {h = c = 
he = 1), where 



1 nm"^ = 197.327 eV, 



(23) 



and thus for a given radius R in nm, the transverse mo 
mentum of a matching Bessel mode in keV is: 

0.37... 



R 



(24) 



Vortex core size. By equating the two terms in (p!6|) 
and expanding the Bessel functions up to O ((/c±r)^), one 
has: 



kl 



jS{k^rc) = J?{k±rc) 



a 1 



2 oh,4 4 



C 



32 



h,'2 '2 7^4 ^4 



16 



(25) 



where a = -^^qz^^^y^pp- • One has to take note that if rc is 
expressed in nm, one has to introduce a conversion factor 
from natural units between the k± in a and the k± in 
(|25]) . The smallest positive root gives an expression for 
small k±: 

; ••• ^ +0(v^)- (26) 
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